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Abstract 

We prove that knowing the lengths of geodesies joining points of the boundary 
of a two-dimensional, compact, simple Riemannian manifold with boundary, we can 
determine uniquely the Riemannian metric up to the natural obstruction. 

1 Introduction and statement of the results 

Let (M,g) be a compact Riemannian manifold with boundary dM. Let d g (x,y) denote the 
geodesic distance between x and y. The inverse problem we address in this paper is whether 
we can determine the Riemannian metric g knowing d g {x,y) for any x G dM, y G dM. 
This problem arose in rigidity questions in Riemannian geometry jM], (Hj, [C^rj . For the 
case in which M is a bounded domain of Euclidean space and the metric is conformal to 
the Euclidean one, this problem is known as the inverse kinematic problem which arose in 
Geophysics and has a long history (see for instance |Rj and the references cited there). 

The metric g cannot be determined from this information alone. We have d^* g = d g for 
any diffeomorphism ip : M — > M that leaves the boundary pointwise fixed, i.e., ip\dM = M 
where Id denotes the identity map and ip*g is the pull-back of the metric g. The natural 
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^Partly supported by NSF and a John Simon Guggenheim Fellowship. 
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question is whether this is the only obstruction to unique identifiability of the metric. It is 
easy to see that this is not the case. Namely one can construct a metric g and find a point 
xo in M so that d g (xo, dM) > sup xye dMdg(x, y). For such a metric, d g is independent of a 
change of g in a neighborhood of xq. The hemisphere of the round sphere is another example. 

Therefore it is necessary to impose some a-priori restrictions on the metric. One such 
restriction is to assume that the Riemannian manifold is simple, i.e., given two points there 
is a unique geodesic joining the points and dM is strictly convex. dM is strictly convex if 
the second fundamental form of the boundary is positive definite in every boundary point. 

R. Michel conjectured in jM] that simple manifolds are boundary distance rigid that is 
d g determines g uniquely up to an isometry which is the identity on the boundary. This 
is known for simple subspaces of Euclidean space (see |Grj ). simple subspaces of an open 
hemisphere in two dimensions (see jM]), simple subspaces of symmetric spaces of constant 
negative curvature [BCG], simple two dimensional spaces of negative curvature (see |Clj or 

10]). 

In this paper we prove that simple two dimensional compact Riemannian manifolds are 
boundary distance rigid. More precisely we show 

Theorem 1.1. Let (M,gi),i = 1,2 be a two dimensional simple Riemannian compact Rie- 
mannian manifold. Assume 

d 9l (x,y) = d 92 (x,y) V(ar, y) G dM x dM 

then there exists a diffeomorphism if) : M — > M, iP\qm = Id, so that 

92 = ip*gi- 

As it has been shown in |Shj . Theorem 1.1 follows from 

Theorem 1.2. Let (M,gi),i = 1,2 be a two dimensional simple Riemannian compact Rie- 
mannian manifold. Assume 

d gi (x,y) = d 92 (x,y) V(x, y) e dM x dM 

and 

gi\dM = 92\dM- 

then there exists a diffeomorphism if) : M — > M, iP\qm = Id, if>'\dM = Id, so that 

92 = ip*gi- 

We will prove Theorem 1.2. The function d g measures the travel times of geodesies joining 
points of the boundary. In the case that both gi and gi are conformal to the Euclidean metric 
e (i.e., (gk)ij = ot^Sy, k = 1,2 with 5y the Kronecker symbol), as mentioned earlier, the 
problem we are considering here is known in seismology as the inverse kinematic problem. In 
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this case, it has been proven by Mukhometov in two dimensions |Mnj that if (M, gi),i = 1,2 
is simple and d gi = d g2 , then g\ = g2- More generally the same method of proof shows 
that if (M,gi),i = 1,2 are simple compact Riemannian manifolds with boundary and they 
are in the same conformal class, i.e. g± = ag2 for a positive function a and d gi = d 92 then 
gi = g2 [Mill]. In this case the diffeomorphism ip must be the identity. For related results 
and generalizations see jB], [EH], P, [GN], [ME]. 

We mention a closely related inverse problems. Suppose we have a Riemannian metric in 
Euclidean space which is the Euclidean metric outside a compact set. The inverse scattering 
problem for metrics is to determine the Riemannian metric by measuring the scattering op- 
erator (see [0])- A similar obstruction occurs in this case with ip equal to the identity outside 
a compact set. It was proven in [G] that from the wave front set of the scattering operator 
one can determine, under some non-trapping assumptions on the metric, the scattering 
relation on the boundary of a large ball. We proceed to define in more detail the scattering 
relation and its relation with the boundary distance function. 

Let v denote the unit-inner normal to dM. We denote by Q (M) — > M the unit-sphere 
bundle over M: 

n(M) = |J n x , Q x = {^e t x (m) ■. \z\ g = 1}. 

xeM 

Q(M) is a (2 dim M— l)-dimensional compact manifold with boundary, which can be written 
as the union dQ (M) = d+Q (M) U9_H (M) 

d ± Q (M) = {(x, e dQ (M) , ± {u (x) , > }. 

The manifold of inner vectors d + Q (M) and outer vectors d_Q(M) intersect at the set of 
tangent vectors 

d n (M) = {( x , o e on (M) , (u (x) ,0 = o}. 

Let (M, g) be an n-dimensional compact manifold with boundary. We say that (M, g) 
is non-trapping if each maximal geodesic is finite. Let (M, g) be non-trapping and the 
boundary dM is strictly convex. Denote by r(x, the length of the geodesic j(x, £, t), t > 0, 
starting at the point x in the direction £ e Q x . These function is smooth on Q(M) \<9of2(M). 
The function r° = r\dn(M) is equal zero on d-Q(M) and is smooth on d + Q(M). Its odd part 
with respect to £ 

r°(x,0 = ^(r°(x,0-r°(x,-0) 

is the smooth function. 

Definition 1.1. Let (M,g) be non-trapping with strictly convex boundary. The scattering 
relation a : dQ (M) — > dQ (M) is defined by 

a(x,0 = (7(x,^,2r°(x,0),7(^e,2r°(x,0)). 
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The scattering relation is a diffeomorphism dfl (M) — > dQ (M) . Notice that at\g + Q(M) '■ 
d+tt(M) -> d-tt(M), a\ 9 _n( M ) ■ d^tt(M) -> d + Vt(M) are diffeomorphisms as well. Ob- 
viously, a is an involution, a 2 = id and <%f2 (M) is the hypersurface of its fixed points, 
a(x,0 = (x,0, (x,0e<9 O(M). 

A natural inverse problem is whether the scattering relation determines the metric g 
up to an isometry which is the identity on the boundary. In the case that (M, g) is a 
simple manifold, and we know the metric at the boundary, knowing the scattering relation 
is equivalent to knowing the boundary distance function ([Mj). We show in this paper that 
if we know the scattering relation we can determine the Dirichlet-to-Neumann (DN) map 
associated to the Laplace-Beltrami operator of the metric. We proceed to define the DN 
map. 

Let (M, g) be a compact Riemannian manifold with boundary. The Laplace-Beltrami 
operator associated to the metric g is given in local coordinates by 



A g u = on M, u 



where (g^) is the inverse of the metric g. Let us consider the Dirichlet problem 

We define the DN map in this case by 

A g (f) = (u,Vu\ 9M ) 

The inverse problem is to recover g from A g . 

In the two dimensional case the Laplace-Beltrami operator is conformally invariant. More 
precisely 

for any function (5, (5 ^ 0. Therefore we have that for n = 2 

for any non-zero (3 satisfying (3\qm = 1- 

Therefore the best that one can do in two dimensions is to show that we can determine 
the conformal class of the metric g up to an isometry which is the identity on the boundary 
That this is the case is a result proven in [LeU] for simple metrics and for general connected 
two dimensional Riemannian manifolds with boundary in [LaU]. 

In this paper we prove: 

Theorem 1.3. Let (M,gi),i = 1,2 be compact, simple two dimensional Riemannian mani- 
folds with boundary. Assume that a gi = a 92 . Then A gi = A g2 . 
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The proof of Theorem 1.2 is reduced then to the proof of Theorem 1.3. In fact from 
Theorem 1.3 and the result of [LaU] we get that we can determine the conformal class of 
the metric up to an isometry which is the identity on the boundary. Now by Mukhometov's 
result we have that the conformal factor must be one proving that the metrics are isometric 
via a diffeomorphism which is the identity at the boundary. In other words d gi = d 92 implies 
that a gi = a g2 . By Theorem 1.3 A gi = A g2 . By the result of [LeU], [LaU], there exists a 
diffeomorphism ip : M — > M, ip\dM — Identity and a function (3 ^ 0, (3\qm = identity such 
that gi = (3ifj*g2- By Mukhometov's theorem (3 = 1 showing that g\ = ip*g2 proving Theorem 
1.2. and Theorem 1.1. 

The proof of Theorem 1.3 consists in showing that from the scattering relation we can 
determine the traces at the boundary of conjugate harmonic functions, which is equivalent 
information to knowing the DN map associated to the Laplace-Beltrami operator. The steps 
to accomplish this are outlined below. It relies in a connection between the Hilbert transform 
and geodesic flow. 

We embed (M, g) into a compact Riemannian manifold (S, g) with no boundary. Let ip t 
be the geodesic flow on Q(S) and 7i = ^v 9 t|t=o be the geodesic vector field. Introduce the 
map ip : fi(M) -> <9_fi(M) defined by 

*p(x, o = ip T {x&)(x, 0, (x, o e n(M). 

The solution of the boundary value problem for the transport equation 

Hu = 0, w|a + n(M) = w 

can be written in the form 

u = = w o a o ip. 
Let u be the solution of the boundary value problem 

7-Lu = -/, u\ d _n( M) = 0, 

which we can write as 

uf(x,£)= J f(<p t (x,Z))dt, (^)Gfl(M). 
o 

In particular 

Ht = -1. 

The trace 

If = u f \ 9+ n( M ) 

is called the geodesic X-ray transform of the function /. By the fundamental theorem 
of calculus we have 

IHf = (f o a - f)\ d+a{M) . (1.1) 
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In what follows we will consider the operator I acting only on functions that do not depend 
on £, unless otherwise indicated. Let L^(<9 + fi(M)) is the real Hilbert space, with scalar 
product given by 

(u,v) L 2 (9+n(M)) = / fiuvdZ, ^ = (f,i/). 

Jd+Cl(M) 

Here the measure dE = d(dM) A dVl x where d(dM) is the induced volume form on the 
boundary by the standard measure on M and 

n 

dQ x = J2{-l) k+1 ^ k d^ A ... A d\ k A ...dC- 

k=l 

As usual the scalar product in L 2 (M) is defined by 

(u,v)— / uv y/detgdx. 

J M 

The operator / is a bounded operator from L 2 (M) into L 2 (d+£l(M)) . The adjoint /* : 
Lj(d+ft(M)) -> L 2 (M) is given by 



I*w{x) = / wj,(x,£)dQ x , 



We will study the solvability of equation I*w = h with smooth right hand side. Let 
w G C°° (d + Q(M)) . Then the function will not be smooth on Q(M) in general. We 
have that G C°°(Q(M) \ d Q(M)). We give below necessary and sufficient conditions for 
smoothness of on Q(M). 

We introduce the operators of even and odd continuation with respect to a: 

A ± w(x, = ± (a*w) (x, 0, Or, e <9_fi (M) . 

The scattering relation preserves the measure |(£, i/)|dE and therefore the operators A± : 
Ll{d + tt{M)) -> (dfi(M)) are bounded, where Lj^ (dfi(M)) is real Hilbert space with 
scalar product 

( w > u )Lf A , | («i(M)) = y N /i =(^,z/). 
9n(M) 

The adjoint of A± is a bounded operator A* ± : (<9fi (M)) -> L^(9 + fi(M)) given by 

A* ± u = (uinoQ)| a+Si(M) . 
Using A I formula (1.1) can be written in the form 

IHf = -A*_f°, f° = f\ d n { M). (1.2) 
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The space (d+Q (M)) is defined by 

(d+Q (M)) = {weC QO (d + Q (M)) : uty e C°° (ft (M))}. 



We have the following characterization of the space of smooth solutions of the transport 
equation 

Lemma 1.1. 

C™{d + Q(M)) = {we C°°{d + Q{M)) : A + w e C°°(9fi(M))}. 

Now we can state the main theorem for solvability for /*. 

Theorem 1.4. Let (M,g) be a simple, compact two dimensional Riemannian manifold with 
boundary. Then the operator I* : C™(d + Q(M)) — > C°°(M) is onto. 

Now we define the Hilbert transform: 



where the integral is understood as a principle value integral. Here _L means a 90° degree 
rotation. In coordinates = where 



The Hilbert transform H transforms even (respectively odd) functions with respect to £ 
to even (respectively odd) ones. If H + (respectively HJ) is the even (respectively odd) part 
of the operator H: 



and u+, U- are the even and odd parts of the function u, then H + u = Hu + , H_u = Hu-. 

We introduce the notation 7ij_ = (£j_,V) = — (£,Vj_), where Vj_ = eV and V is the 
covariant derivative with respect to the metric g. The following commutator formula for the 
geodesic vector field and the Hilbert transform is very important in our approach. 

Theorem 1.5. Let (M,g) be a two dimensional Riemannian manifold. For any smooth 
function u on Q(M) we have the identity 




(1.3) 



e = 





[H, H]u = n±u + (Wju)o 



(1.4) 



where 




is the average value. 
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Now we can prove the Theorem 1.3. 

Separating the odd and even parts with respect to £ in (1.4) we obtain the identities: 

H + Tiu — HH-u = (Tt±u)o, HJHu — HH + u = H±uo. 

Let (M,g) be a non-trapping strictly convex manifold. Take u = w^,w G C£°(<9+(Q)). 
Then 

2nTCH + w,p = —H±I*w 
and using the formula (1.2) we conclude 

2nA*_H + A + w = IH ± I*w, (1.5) 

since w^\ dn ^ M ) = A + w. 

Let (h, h*) be a pair of conjugate harmonic functions on M, 

V/i = Vj./i*, V/i* = -Vj./i. 

Notice, that <5V = A is the Laplace-Beltrami operator and 5V± = 0. Let I*w = h. Since 
IH±h = ITih^ = —A*_h®, where h® = hJg M , we obtain from (1.5) 

2ttA*_H + A + w = -A*_h°. (1.6) 

The following theorem gives the key to obtain the DN map from the scattering relation. 

Theorem 1.6. Let M be a 2- dimensional simple manifold. Let w G C^(d + Q(M)) and /i* 
is harmonic continuation of function h®. Then the equation (1.6) holds iff the functions 
h = I*w and h* are conjugate harmonic functions. 

Proof. The necessity has already been established. Using (1.2) and (1.5) the equality (1.6) 
can be written in the form 

mji = IHq, 

where q is an arbitrary smooth continuation onto M of the function h® and h = I*w. Thus, 
the ray transform of the vector field Vg + V±h equal 0. Consequently, this field is potential 
([An]), that is, Vg + V±h = Vp and p\ou = 0. Then the functions h and h* = q — p 
are conjugate harmonic functions and hJgM = h®. We have finished the proof of the main 
theorem. 

□ 

In summary we have the following procedure to obtain the DN map from the scat- 
tering relation. For an arbitrary given smooth function h® on dM we find a solution 
w G C^(d + Q(M)) of the equation (1.6). Then the functions h° = 2tt(A + w) (notice, that 
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2ir(A+w)o = I*w\qm) and h® are the traces of conjugate harmonic functions. It gives the 
map 

which is the DN map. This proves Theorem 1.3. For the definition of V" see section 2. 

A brief outline of the paper is as follows. In section 2 we collect some facts and definition 
we need later. In section 3 we study the solvability of I*w = h on Sobolev spaces and prove 
Theorem 1.4. In section 4 we make a detailed study of the scattering relation and prove 
Lemma 1.1. In section 5 we prove Theorem 1.5. 

2 Preliminaries and Notation 

Here we will give some definitions and formulas, that we need in what follows. For further 
references see jEj, jJ], jKj, |Shj . Let it : T(M) M be the tangent bundle over an n- 
dimensional Riemannian manifold (M,g). We will denote points of the manifold T(M) by 
pairs (x, £). The connection map K : T(T(M)) — > T(M) is defined by its local representation 

K(x, £, y, V ) = (x, V + T(x)(y, £)), (T(x)(y, £))* = V) k {x)y^\ i = 1, n, 

where T l - k are the Christoffel symbols of the metric g 

p 1 n ii( d 9ji | d 9ki dg jk 
3k 2 y { dx k dxi dx l 
The linear map K(x,£) = K\i x n : Tt X) n(T(M)) — > T X M defines the horizontal subspace 
H( Xi £) = KerK(x,^). It can be identified with the tangent space T X (M) using the operator 

The vertical space V< x £) = KerTr'(x,^) can also identified with T X (M) by using the operator 

J^ = {K{x,t)\v M )- 1 ■■ T X {M) -> V M . 

The tangent space T( X £)(T(M)) is the direct sum of the horizontal and vertical subspaces, 
T( X: ^)(T(M)) = Hfag) © V(x,£)- An arbitrary vector X G T^ X ^(T(M)) can be uniquely 
decomposed in the form 

X = J( x £)Xh + J( x £)Xv, 

where 

X h = ir'(x,OX, X v = K{x,i)X. 

We will call Xh,X v the horizontal and vertical components of the vector X and use the 
notation X = (Xh,X v ). If in local coordinates X = (X 1 , ...,X 2n ) then X^,X V is given by 

X\ = X\ X l v = X l+n + T) k (x)X^ k , i = l,...,n. 
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The horizontal and vertical components Y h ,Y v G T*(M) of a covector Y G T* x ^(T(M)), Y = 
(Y h , Y v ) are given by 

(X,Y) = (X h ,Y h ) + (X v ,Y v ). 

In local coordinates we have 

Yj 1 — Yj — T l jk (x)Y i+n £ k , Y? = Yj +n , j = 1, ...,n. 

Let iV be a smooth manifold and / : T(M) — > iV a smooth map. Then the derivative 
f'(x,£) : T( Xj £)(T(M)) — > Tf( X £)(N) defines the horizontal Vhf(x,Q and vertical V v f(x,£) 
derivatives: 

V fc /(x,0 = o 4, : T X (M) - T /(Xi€) (JV), 

V w /(x,0 = o J ( ^ } : T X (M) -> T /(Xi0 (JV). 

We have then that 

/'Or, OX = (V h /(x, 0, X,) + (V„/(x, 0, X.) (2.1) 

In local coordinates 

v w / (a) (*,0 = ( A _ r^)^)/^^), 

V Bi / (a) (*>0 = A/Mfofl, a = l,...,dimiV. 

Let e : £ — > T(M) be a smooth embedding and p : T(M) — > £ denotes a normal 
projection, i.e. pe = id. For a point (a;,f) G e(E), Kerp'(x,f) = JV(E) f| T (x ^(T(M)), 
where X(£) is the normal bundle of E. 

Definition 2.1. TTie horizontal and vertical derivatives Vj[, V|J on E are £/ie tangent com- 
ponents ofVh and V v with respect to E: 

vi/ = V fc (/op)| E , Vj/ = V„(/op)| E . 

This definition is obviously independent of the choice of the normal projection. From the 
identity / = (fop) o e and (2.1) we have for any vector X G T( Xj £)(T(£)), by identifying X 
with its embedding e'(x,£)X, that 

f'(x, OX = (/ o p)'(x, OX = (V h (/ o p)(x, 0, X,) + (V„(/ o p)(x, o, x,). 

Therefore 

/'(x,ox = (vlf(x,t),x h ) + (vil/(x,o,x,). 



10 



Let us give an equivalent definition. Denote by i(x,£) the isomorphism 

% (x, : T (x£) (T (M)) - T X (M) x T x (M) , defined by i (x, £) X = (X,, X„) 

Let is be its restriction on T( x ,£) (E) , is (x, £) = i (x, £) |r. s) (s)- Denote by T x ^ (M) = 
prj Ranij] (x,£) C T X (M) x T s (M) , i = 1,2, the projections to the first and second com- 
ponent of the range of is- Then 

Vj[/(s,0 = V fc /(x,0 | T (D (M) , Vj/(x,0 = V„/(a;,0 | T ( 2)(M) , 

where / : T (M) — > X is an arbitrary smooth continuation of /. 
Example 1 

Let us consider the special case of £ = Q(M) the unit sphere bundle. Then p(x,£) = 
(x, is a normal projection. For simplicity of notation in the case of the manifold Q(M) 

we will use notation d instead of Vll. Then 

9 = (^-£(£>^))Iki=i- 

Clearly, (£, <9) = 0. Notice, that the horizontal derivative is tangent to the submanifold Q(M) 
since V/J£| = 0. Because of this we keep the original notation for horizontal derivative on 
Vl(M) and moreover simplify this notation to just V. Thus, in the case of the manifold f2(M) 
we will use notation V and d for the horizontal and vertical derivatives. We will also use V 
for horizontal derivative on T(M). Notice relation 

7W)(Q(M)) = {X G T ( ^)(T(M)) : (£, JQ = 0}. 

Example 2 

Let T be a smooth hypersurface in M with normal v and E its lift to Q(M), E = 7rf 1 (r) = 
{(x, £) :xGr,(G ^(M)}, where 7Ti = 7r|n(M)- Then 

^(^(r)) = {X G T ( ^)(T(M)) : (i/(x),X fc ) = 0, (Z,X V ) = 0} (2.2) 

It means that TP (M) = {rj G T x (M) : (rj, 1/ (a;)) = 0}, Tj 2) (M) = {77 G T x (M) : (£,77) = 
0}. We have the same vertical derivative as in the case of Q (M) . The horizontal derivative 
for this case will be denoted by V", so that 

v» = (v-i/(i/,v))L r i (r) . 

Clearly, (u, V 11 ) = 0. 

We now state a similar definition of vertical and horizontal derivatives on a submanifold 
for semibasic tensor fields. We recommend chapter 3 of [Sh] for more details. 
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Let TJ(M) denote the bundle of tensor fields of degree (r,s) on M. A section of this 
bundle is called a tensor field of degree (r,s). Let ix r s : TJ(M) — > M be the projection. 
A fiber map u : T(M) — > T^(TM), i.e. 7r£ o -u = 7r is called a semibasic tensor field of 
degree (r, s) on the manifold T(M). Denote by £ the semibasic vector field given by the 
identity map T(M) — > T(M). An arbitrary tensor field w of degree (r, s) on the manifold 
M, i.e. section -u : M — > TJ(M) defines by the formula w o 7r a semibasic tensor field (since 
7r£ o (-u o 7r) = {-K r s o -u) o Ti = id o 7r = 7t). The map -u — * -u o 7r identifies tensor fields on M 
and ^-constant semibasic tensor fields on T(M). A semibasic tensor field on the submanifold 
E C T(M) is defined as the fiber map u : E — > TJ(M), i.e. 7r£ o -u = 7r| s . If w is semibasic 
tensor field on T(M), then its restriction w| s is a semibasic tensor field on E. Using the 
metric g we can identify the bundle TJ(M) with T r+S (M) and the bundle T r ° +S (M) with 
T r * +S (M). 

A semibasic tensor field u : T(M) — > T^(TM) is naturally identified with the polynomial 
of degree m in rj, u (m) : T(M) x T m (M) -> i?, w m (x, £, 77(1), ?? (m )) = w^...^ (x, ?7(i)-»7(m)' 
where the fibers of the vector bundle T(M) x T m (M) = T m+1 (M) are given by 

T x n+1 (M) = T x (M)x . . .xT x (M) . 

^ V ^ 

The horizontal derivative V/iW( m ) and vertical derivative V„M( m ) are the fiber maps 
T(M) x T m (M) -> T*(M) = T°(M), and therefore are semibasic tensor fields on T(M) x 
T m (M) of degree 1. The corresponding polynomials (V/ l «( m ))(i), (V^(m))(i) on T(M) x 
T m+1 (M) define semibasic tensor fields Vw, V^w of degree (m + 1), which we call the hori- 
zontal and vertical derivatives of u. Thus, by definition 

(V«)( m +1) = (V h U( m ))(l), (V ? «)( m+ l) = ( V^(m)) (1) . 

Straightforward calculations give the representation in local coordinates: 

\ V a )h...i m+1 - ^im+i a h-im i i (m+1) fc? Qtj ' l V ^J»l..-Wl - 0£i m+ i ' 

where V denotes the usual covariant derivative on the manifold (M,g). Notice, that for 
^-constant tensor fields, Vw = Vw and since we identify ^-constant semibasic tensor fields 
with tensor fields on M, we will use one notation V for covariant and horizontal derivatives. 

As in the case of maps, we define tangent derivatives of semibasic tensor fields on the 
submanifold £ C T(M): 

V ll w = V(uop)| E , Vjw = V e (uop)| s . 

In the case £ = f2(M) we keep the notation V instead of V" and use d instead vjj. We 
mention the following formulas, [Sh] : 

Vg = 0, V£ = 0, 0^-8]- Ctr (2.3) 
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[V,0] =0, ();. <)j 0<h-C/h. 



(2.4) 



[V t ,V J }u = -R p qiJ d p u, 



(2.5) 



where i? is the curvature tensor. In the last formula u is a scalar. In the case that S = 7if 1 (r) 
we will use the same notations V", d as in the case of maps. 

3 The geodesic X-ray transform 

In this section we study the solvability of the equation I*w = h and prove Theorem 1.4. 

Lemma 3.1. Let V be an open set of a Riemannian manifold (M,g). We can define the ray 
transform as before. Then the normal operator 1*1 is an elliptic pseudodifferential operator 
of order —1 on V with principal symbol c n where c n is a constant. 

Proof. It is easy to see, that 



(/*//) (a;) = dn x / / ( 7 (x, C, t)) dt = 2 / dn x / / ( 7 (x, £, t)) dt. (3.1) 



Before we continue we make a remark concerning notation. We have used up to know 
the notation 7(2, £,t) for a geodesic. But it is known jj] , that a geodesic depends smoothly 
on the point x and vector ^ 6 T X (M). Therefore in what follows we will also use sometimes 
the notation 7(2;, £i) for a geodesic. Since the manifold M is simple and any small enough 
neighborhood U (in (S,g)) is also simple (an open domain is simple if its closure is simple). 
For any point x G U there is an open domain C T x (U) such that exponential map 
exp x : — > U, exp x r] = 7(2, rj) is a diffeomorphism onto U. Let D x , x G M be the inverse 
image of M, then exp x (D x ) = M and exp x \o x : D x —>■ M is a diffeomorphism. 

Now we change variables in (3.1), y = j(x,£t). Then t = d g (x,y) and 






M 



where 



K{x,y) = 2 



det (exp x 1 ) / (x, y) y/ det g (x) 
d^ 1 (x,y) 



Notice, that since 



7(2,77) = x + r] + 0(\r]\ 2 ), 



(3.2) 
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it follows, that the Jacobian matrix of the exponential map is 1 at 0, and then det(exp x 11 )(x, x) = 
1/det (exp x ) (x, 0) = 1. ^From (3.2) we also conclude that 

d 2 (x, y) = (x, y) (x - yf (x - y) j , G tj (x, x) = 9ij (x) , G tJ G C°° (M x M) 

Therefore the kernel of 7*7 can be written in the form 

2 det (exp' 1 )' (x, y) yj det g i 



K(x,y) ~ ' " 



(Gij{x,y) (x - y) 1 {x - yf^j 



(n-l)/2- 



Thus the kernel K has at the diagonal x = y a singularity of type \x — y\ n+1 . The 
kernel 

K (xv)- 2y/detg(x) 

^0 \F,V) ~ , _ \ (n-l)/2 



gtj (x) (x - y) 1 {x - y)^ 



has the same singularity. Clearly, the difference K — K has a singularity of type \x — y\~ n+2 ■ 
Therefore the principal symbols of both operators coincide. The principal symbol of the 
integral operator, corresponding to the kernel Kq coincide with its full symbol and is easily 
calculated. As a result 



a (1*1) (x, = 2 v/det^M / --^^—-dy = c n ^ . 



e -i(y,0 
{g l3 {x)y^t- l)/21 

□ 



Let r M denotes the restriction from S onto M. 

Theorem 3.1. Let U be a simple neighborhood of the simple manifold M . Then for any 
function h G H s (M) , s > there exists function f G -£P _1 (U) , r M I*If = h. 

Proof. Let (M, g) be simple and embedded into a compact Riemannian manifold (S, g) 
without boundary, of the same dimension. Choose a finite atlas of S, which consist of simple 
open sets Uk with coordinate maps n k : U k — > R n . Let {<Pk} be the subordinated partition of 
unity: ipk > 0, suppifk C Uk, YliVk — 1- We assume without loss of generality that M C U± 
and </?i]m = 1- We consider the operators Ik, II for the domain Uk, and the pseudodifferential 
operator on (S, g) 

k 

Every operator 7^7^ : (Uk) — > C°° (C4) is an elliptic pseudodifferential operator of order 
— 1 with principal symbol c n , £ G T (Uk) ■ Then P is an elliptic pseudodifferential 
operator with principal symbol c n , £ G T (£) and, therefore, is a Fredholm operator 
from H S (S) into H S+1 (S). We have that Ker P has finite dimension, Ran P is closed and 
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has finite co dimension. Notice, that P* = P (more precisely if P s = P : H s (S) — > H s+1 (S) , 
then (P s )* = P_ a _x). 

For arbitrary s > the operator tm '■ H s (S) — > P s (M) is bounded and tm{H s (S)) = 
H s (M) . Then the range of r M P : H s (S) -> P s+1 (M) , s > -1 is closed. 

Since M is only covered by U\ and </?i|m = 1 we have that vmPJ = tmP\Ii (fluj- Thus, 
the range of the operator r M I{Ii : H s (Ui) — > P s+1 (M) , s > — 1 is closed. Now to prove 
the solvability of the equation 

r M /*P/ = h G H s+1 (M), s> -1 

in H s (U ± ) is sufficient to show that the kernel of the adjoint {r M I{h)* ■ (P (s+1) (M))* -> 
(P s (Pi))* is zero. 

Let (,) M and (,) be dualities between H S (M) and (P s )* (M) or H S (S) and H- S (S) 
respectively. The dual space (H s (M))* , s > can be identified with the subspace of 
H~ s (S):' 

(H s (M))* = H- s (M) = {ue H- s (S) : suppu C M) . 
For any / G H s (U{) , if G P~( 1+s ) (M) we have 

(r M Ilhf,u) M = (P s f,u) = (/,P_ s _iu), 

where / is an arbitrary continuation of / on the manifold S. On the other hand 

{r M I{hf,u) M = (/, {r M Ilh)*u) M . 

Since / is arbitrary, then equality (/, P_ s _iu) = (/, (rj^-P^P)* u)m implies (rjyP^P)* = 
Tux P -s-\ = ruJih- 

Because of ellipticity the equality rjj 1 Pu = implies smoothness w^, and then w G 
P^ 1 (M) implies of w G C °°(Pi). Since r^Pu = J*Pm, then J*Pm = =>• iiu = =>• 
« = 0. 

□ 

Now we are ready to prove Theorem 1.4. 

Proof. Let /, p be the geodesic X-ray transforms on M and P x respectively. From Theorem 
3.1 it follows that for any h G C°°(M) there exists / G C°°(P), such that r M I{hf = h. 
Then G C°°(Q(7V)) where TV denotes the interior of N. Let u> = 2u+\d + n(M), where u+ is 
the even part with respect to £. Then it easy to see, that it;^, = 2tt+|^(M) and I*w = h. The 
function w G C~(d+Q(M)) since G C°°(fi(M)). □ 
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4 Scattering relation and folds 



In this section we prove Lemma 1.1. 

As indicated before we embed (M, g) into a compact manifold (S, g) with no boundary. 
Let (N,g) be an arbitrary neighborhood in (S,g) of the manifold (M,g), such that any 
geodesic 7(2;, £,t), (x,£) G Q(N) intersects the boundary dN transversally. Then the length 
of the geodesic ray r is a smooth function on Q(N) and the map ifj : dCl(M) — > d+Q(N), 
defined by 

is smooth as well. Moreover it turns out is a fold map with fold d £l(M). This fact will 
be proved in the next Theorem. Once this is proven Lemma 1.1 follows from [H] Theorem 
C.4.4. From the assumption A + w G C°°(dQ(M)) we deduce the existence of a smooth 
function v on a neighborhood of the range <p(dQ(M)) such that w — vocf). Consider function 
= w o a o ip. Change notation ip to tpM, keeping w^. Denote by ipN the map, analogical 

tO lf) M , 

i>N (X, = ip T ( x £) (x, , (x, g Q (N) . 
Then w$ = v o (ft o a o ip M . It easy to see, that <fi o a o ip M = iPn\q(m)- Since the map ip^ 
is smooth on Vt (M) , then E C°° (SI (M)), i.e. w e C™(d+VL(M)). Thus Lemma 1.1 is 
proven once we show that is a fold. 

Theorem 4.1. Let (M,g) be a strictly convex, non-trapping manifold and N an arbitrary 
neighborhood of M, such, that any geodesic 7(2;, £, t), (x, £) G Q(N) intersects the boundary 
dN transversally. Then the map <fi, defined by (4-1) is a fold with fold doQ(M). 

First we recall the definition of a Whitney fold. 

Definition 4.1. Let M, N be C°° manifolds of the same dimension and let f : M — >■ iV be 

a C°° map with f(m) = n. The function f is a Whitney fold (with fold L) at m if f drops 
rank by one simply at m, so that {2;; df(x) is singular } is a smooth hypersurface nearm and 
ker(df{m)) is transverse to T m L. 

We rewrite the definition of a fold below in a form that we will use. 

Let 7 : (a, b) — > M, G (a, b) be an arbitrary curve in M. Then / induces the curve 
71 : (a, b) — * N in the manifold N, 71 = / o 7. The derivative f'(m) of the map / at the 
point m maps the tangent vector £ = 7(0) G T m (M) into the tangent vector r\ = 71 (0) = 
f'(m)£ G Tf( m (N). The (acceleration) vectors 7(0), 71 (0) belong to the spaces T^(T m (M)) 
and T v (Tf( m )(N)) respectively, which are tangent to the vector spaces T m (M) and Tf( m )(N) 
and may be identified with them. Fix a point m G M and nonzero vector £ G Kerf'(m). 
Let us consider the map L^ : T m — > 7/( m ), ^(7(0)) = 71 (0). A calculation in local 
coordinates gives that the range of the map L^ is an affine subspaces in Tf^(N), RanL^ = 
a(£) + Ranf'(m), where a(£) is a vector, depending on £. Moreover the function 

Hessf(m)(C,Y) = (L^(X),Y),Y G (.Ranf (m)) x 
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does not depend on X. Here (, ) corresponds to the duality of the spaces T n (N) and T*(N). 
It follows immediately from the that the definition of a fold is equivalent to 

Definition 4.2. The map f is a fold at the point m if 

(i) dimKer f'{m) = 1. 

(ii) Hessf(m)(C,Y)^0. 

Remark 4.1. Condition (i) implies that dimCoker f'(m) = 1 and Ran is an hy- 
perplane. and condition 4.2(H) implies that the hyperplane RanL^ does not pass trough the 
origin. 

The following proposition is immediate. 

Proposition 4.1. Let f : M — > N be a smooth map, dimM = dimiV. Let e : M — > M be 

an hypersurface and the map foe: M — > / (M ) is non-singular. Then dim Ker f (x) < 1 
for any point x G M ■ 

Now we start proving Theorem 4.1. 

Proof. We first prove Condition 4.2(i). The smooth map 4>\d n(M) '■ d Q(M) — >■ (j)(d fl(M)) 
is bijective and therefore dimKer 0'(x o ,Co) < 1 for any point (xo,£o) e d Q(M). We prove 
that the vector X G T( XOj £ )(cft)(M)), X oft = Co 5 A" t, — 0, where X 0/ j,X t, are its horizontal 
and vertical components (see section 2) belongs to Ker (j)'(xo,£o). Represent = ip^ ° e, 
where e : dQ(M) — > fi(iV) is the embedding. Then (f>'(x,£)X = i[)' N (x,£)X (we identify here 
the vector X with its embedding e'(x,C)X e 7W)(<9fi( M )))- Then e ifer 0'(^o, Co) iff 
X G T (x0i?o) (9fi(M))f|^er^(x ,^ )- The fact' that X G T (ie0i&) (dfi (M)) follows from 
equality (^(z ),Co) = (see (2.2)). The equality ^atO^o, Co)^o = (V^a^o, Co), Co) = is 
obvious, since for any point (x,£) G Q(N) the point iPn(x,^) is left fixed by geodesic flow 
and (V^at^, C), C) = Hi) n = 0. Thus, the vector X = (Co,0) G Ker 0'(xo,£o) an d 
the first condition in the definition of a fold has been verified. 

Let us prove the second condition in the definition of a fold by contradiction. Assume 
that Hess <j>(x , £o)(X , Y°) = 0, where Y° G Ker(<j/(x , Co))*- We will show then Y° = 0. 
Denote (y, rf) = 4>(x, £). Then in local coordinates we have for the covector Y° = (ly , Y^): 

Yyly\x , Co) + Y° +l Vl V \x , Co) = 0, (4.2) 

Yfdktfixo, Co) + Y^d^ixo, Co) = 0. (4.3) 

We write down these equations using the horizontal and vertical parts of the covector Y° (see 
section 2) and the Jacobi fields (x ,£ ,t) , B^) ( x o,£o,t) , k — 1, ...,n on the geodesic 
7(x , Co, t) : 

A*0 = V fe7 + 7V][r, 5 fe = <9 fc 7 + 7<9 fe r. 
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Then the equations (4. 2), (4. 3) can be written in the form 

[(A (k) X) + (DA {k) X)]\t=r = 0, 
[(B (k) X) + (DB ( k)X)]\t=T = 0, 

where D means covariant derivative along geodesic. Without lost of generality one can take 
Yy° = Z c \ t=T , Yfi = —DZ c \ t=T , where Z c is some Jacobi field on j(xq, £o,t). It is known that 
for any Jacobi fields X, Y, the Wronskian {X, Y} = (DX, Y) - (X, DY) is constant. Then 
we get 

{Z c ,A {k) }\ t=0 = 0, {Z c ,B (k) }\ t=0 = 0. 
The initial data is given by 

A J {k) \ t =o = 8{ - u l u ok + CoVjJr, DA (k) \ t=0 = 0, 

5(fc)|t=o = £o<9fcT, (-D-B(fc))- 7 |t=o = - £o£ofc, 
where u = v{xq). We have 

Z c \ t=0 = a£ + cdr, DZ c \ t=0 = bu - cV ll r, 

where a, b, c are some constants, depending on the point (xq, £q). From these general solutions 
we find one satisfying the condition that Y° e T^ x ^(dft(N)), or the conditions (2.2), 

X, Mv)) = -(Dz c \ t=T , My)) = o, (4.4) 

X:V) = (Z C \t=r:V)=0: (4-5) 

where is the normal to ON . Consider at first the condition (4.5). Using the Jacobi field 
J = (r — t)i it can be rewritten in the form 

{Z c ,J}\ t=T = {Z c ,J}\ t=0 = 0, 

and, therefore, a = —cr. Thus, 

Z c \ t=0 = c(dr - Cor), DZ% =0 = bu - cV"r. 

Now consider our hypothesis that Hess 0( XO ^ O )(X O , Y°) = 0. A straightforward calculation 
shows that in local coordinates we have 

Hess <f> M) (X X) = -B(x ,Co){Z,P}\t= T = -B(x ,^){Z c , P}\ t=0 , 

where B = — ^^ Vf^j is the second fundamental form of the boundary DM and 

P = {y, V 7 ) + iy, Vr)7 
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is a Jacobi field. Since B > 0, our hypothesis gives that {Z, P}\ t=0 = 0. Because of the 
initial data P(xo, £o, 0) = vq + (z^o, Vr)£, DP(x , £ , 0) = we obtain 6 = — c(i/ , Vr) and 

= c(dr - Cor), L>Z c | t=0 = -cVr. 

Now we prove that for any point (x, £) G fi(iV) the vector field 

Z (x, £, = (9r - £r) (x, £)) , t e [-r (x, -£) , r (x, 0] (4.6) 

on the geodesic 7 (x, £, t) is a Jacobi field and 

DZ(x,£,t) = (Vr)(^(x,0). (4.7) 

Applying to the equation Hr = —1 the operator rid and using (2.3)-(2.5) we get TC 2 X + 
RX = 0, (AX - )* (x, = i2j fcJ (x) where the semibasic vector field X = dr - f r. The 

operators H and D are related by: 

(HX)o<p t = D(Xo<p t ), (4.8) 

where X is an arbitrary semibasic vector field. Then (4.6) is a Jacobi field. Since from 
Hr = -1 it follows H(dr - £r) = - Vr, we obtain (4.7) from (4.8). 
Thus, Z c (x ,^ ,t) = cZ(x ,!;o,t) and condition (4.4) gives 

= (DZ\ t=T , v N {y)) = -c(Vr(y, V ), v N {y)). (4.9) 

Now we show that 

VT\ d _ n(N) = (4.10) 

Let /i be a smooth function on N, such that h\dM = 0, V/i|m 7^ 0. Then since /i(x, £, r(x, £)) = 
0, (x,£) G Sl(iV), we have that 

v?(y) V7( fc) (a:, r(x, 0) + K (y) , 77) V fe r(x, = 0, 

where (y, 77) = </? T ( x ^)(x, £) G d-Cl(N). Change here (x, £) to </? t (x, £), t G [0, r(x, £)] . Then 
the point (y, 77) does not change. Using the initial data V7^(x, £, 0) = 5 l k we get in the 
limit isn(v) + (^jv(y), ^)Vr(y, 77) = 0. Thus, (4.10) is proved and we conclude from (4.9) that 
c = 0. 

Thus, Hess (;EO)Co) (Xo, Y°) = iff Ker (0'(x o ,£ o ))* = {0}. We have finished the proof 
of the Theorem. 

□ 
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5 The Hilbert transform and geodesic flow 



In this section we prove Theorem 1.5 in the introduction. 
Let H be the Hilbert transform as defined in (1.3) 

We have that H is a unitary operator in the space Lq(Q x ) = {u G L 2 (Q X ) : uq = 0}, 

(u,v) = (Hu,Hv), Vu, v G Lq(Q x ), 

H 2 (u) = -u, Vu G L 2 (Q X ). 

Clearly, all these properties remain the same if we change Q x to f2(M). 

In order to prove Theorem 1.4 we need the following commutator formula which is valid 
for Riemannian manifolds of any dimension 

Lemma 5.1. Let u be a smooth function on the manifold Q 2 (M) = \J x£M Q x , Q x = 
{{x,£,rj) : £,r] G Q x }. Then 

V J u(x, £, v)dn x (77) = J V^u(x, £, 77)^ (77) , (5.1) 
where V^ 2 - 1 under the integral sign in (5.1) denotes the horizontal derivative on VL 2 (M), 

Notice that the horizontal tangential derivative can be defined on T(M) x T(M) in a similar 
fashion to the case of T(M) in section 2. 

Proof. Let G Cq° (R + ) be arbitrary function. We define the function v on T 2 (M) by 

v(x,£,ri) =<p(\ri\)u(x,Z/\Z\,ri/\ri\) 

Let us consider the integral 

S(x,0= f v(x^ : rj)dT x (rj). 

T X (M) 

Identifying T X (M) with R n we have 



S(x,£) = J v(x,£,ri)y/detg(x)dri. 
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Then 



ds ^ ,os r,dv 



™ >t 9v . n . f d In J det g (x) n , 



R" R n 

Since <91n ^det g/dx j = Y k - k we rewrite the last integral in the form 
Then 



R" 

Since 



then after changing to spherical coordinates we obtain 

00 

VS(x, = J V {t) t n ~ l dt J Vu(x, £, v)dn x (77) . (5.2) 



Si, 

Now S in spherical coordinates is given by 



00 

S(x, = y v? (0 r- 1 ^ y £, 77)d«x (v) (5.3) 



fix 

We conclude (5.1) using (5. 2), (5. 3). 

Now we prove Theorem 1.4. 
Proof. A straightforward calculation gives that 

y 1 + «■"> =<) 

and therefore we have 



VHu(x, = ^ / l±(il^) V «(x, 77)^ (77) . 



f2x 
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□ 



For any pair vectors £, r] G £l x we have 



7/ 



Then 



7/ 



i + fog) _ AZ,v) + (Z,v) 2 



Thus 



= HHu + H±u + (H±u) 



and Theorem 1.4 is proved 



□ 
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